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Abstract 

Cd 

We determine the weak mixing angle in the MS-scheme, sin 9w{u), at energy scales 

C^- \ u relevant for present and future low energy electroweak measurements. We relate the 

renormalization group evolution of sin 2 6w( a ) to the corresponding evolution of a{u) and 
. include higher-order terms in a s and a that had not been treated in previous analyses. We 

also up-date the analysis of non-perturbative, hadronic contributions and argue that the 
associated uncertainty is small compared to anticipated experimental errors. The resulting 
value of the low-energy MS weak mixing angle is sin 2 6w{fy = 0.23867 ± 0.00016. 
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§ 1 Introduction 
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The weak mixing angle is one of the fundamental parameters of the electroweak Standard Model 
(D ' (SM). It can be defined through the relation, 

L" 1 /2 

m-k sin 2 W = _S , (1) 

^ ■ 9+9 

where g and g' are the SU (2) L and U (l)y gauge couplings, respectively. Its value is not predicted 
and needs to be extracted from parity violating neutral current experiments, where by far the 
most precise results were obtained at the Z factories LEP 1 and SLC Electroweak symmetry 
breaking provides masses for the W and Z bosons proportional to their gauge interactions. 
Therefore, one has the additional relation, 

sin 2 ^ = l-^|, (2) 

and the gauge boson mass ratio provides independent precise information on sin 2 9w- Extracting 
the fine structure constant, 

e 2 g 2 sin 2 9 w 
a = ~r~ = a > ( 3 ) 

4:71 All 

from the quantum Hall effect or the anomalous magnetic moment of the electron, then fixes both 
gauge couplings. Eqs. p]) and fl3J) are valid at the tree level and modified by radiative corrections. 
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As a result, the precise numerical value of sin 2 6\y depends on the renormalization scheme and 
scale chosen. For example, the on-shell renormalization scheme promotes Eq. (|2"1) to a definition 
of sin 2 9w = sin 2 6^ n_she11 to all orders in perturbation theory. This definition has the advantage 
of being directly related to a physical observable — but only to one-loop order. Since the gauge 
bosons are unstable particles their masses become ambiguous starting at two-loop precision. 

As an alternative, one can define flavor- dependent "effective" mixing angles appearing in the 
Z vector coupling^, 

v f = T f -2Q f sm 2 9f, (4) 

where Qf and Tf are the fermion charge and third component of isospin, respectively. Gauge 
boson self-energy and Zf /-vertex corrections are absorbed into scheme dependent form factors, 
Pf and Kf that are equal to unity at tree-level. The pf are corrections to the overall coupling 
strengths and the Kf are defined by, 

sin 2 Of = K f sin 2 9™~ shen = k f sin 2 9 W (M Z ), (5) 

where the caret marks quantities in the modified minimal subtraction (MS) scheme [2H1]. This 
effective sin 2 9f S is a useful definition as long as electroweak box contributions can be neglected; 
since these do not resonate this condition is clearly satisfied at LEP 1 and SLC However, with 
the greater precision that could be achieved with the GigaZ option of TESLA, such boxes could 
become non-negligible. Thus, it is not easy to construct a definition that can be equated to 
a physical observable to all orders. Neither is this of practical relevance: as long as it is well- 
defined, sin 2 9w can be looked upon as a mere bookkeeping device and means to compare various 
experimental results. 

What is of practical importance, however, is the numerical value of the mixing angle used 
in computing a given observable in a specified renormalization scheme. Generally speaking, one 

2 

expects a one-loop radiatively corrected result to be valid up to small corrections of 0{ ^ 2 si "4 e ) ~ 
10 -4 . On the other hand, for leptons we have, 

(ki - l) 2 « 1.5 x 10~ 3 , (6) 

which is not much smaller than a typical one-loop contribution. The reason is that K£ contains 
top-quark mass enhancements factors, m 2 /M^, that spoil the expected behavior of the pertur- 
bation series. Fortunately, such enhancement factors can be avoided by a judicious choice of 
definition of the weak mixing angle (or renormalization scheme), rendering the truncation error 
small. The MS-definition considered in this article has this property, except that small In m 2 /M^ 
corrections cannot be decoupled simultaneously from all observables [5]IS] . 

For processes off the Z-pole, enhancement factors of similar magnitude as those entering 
Eq. can arise from large logarithms lnM|/m 2 , where mj is some fermion mass. These 
can occur even within a specifically chosen renormalization scheme. Typically, such logarithms 

1 We normalize Vf without an additional factor inversely proportional to sin 9\y cos 6w ■ The scale in this factor 
needs to be chosen of electroweak size, \x = Mz, even for low energy processes. This is automatically achieved 
by normalizing neutral-current amplitudes using p Gf, where Gf is the Fermi constant and p the low energy 
neutral-current p parameter jlj (which is free of fcrmionic mass singularities). Thus, this factor does not affect 
our discussion. 
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are artifacts of using a value for sin 2 9\y obtained at the Z scale (where as noted above, it is 
measured very precisely) in theoretical expressions for very low (or high) energy observables. In 
some renormalization schemes, the weak mixing angle depends explicitly on a renormalization 
scale parameter, which could be the t'Hooft scale, /i, appearing in dimensional regularization, 
or the momentum transfer q 2 . For others there is no explicit scale parameter, and the scale 
dependence is of indirect nature. In either case, the aforementioned logarithms are a potential 
hazard and should be avoided or re-summed if possible. 

The main goal of this paper is to present an analysis of the weak mixing angle in the MS- 
scheme, sin 2 9^(0), relevant for observables measured at (almost) vanishing momentum transfer. 
These observables include the nuclear weak charge obtained in the well-known cesium atomic 
parity violation measurement performed by the Boulder group [TJ; the parity violating M0ller 
asymmetry at SLAC [5] ; and the up-coming measurement of the weak charge of the proton at Jef- 
ferson Lab [H] . Precision measurements of the parity violating deep inelastic eD asymmetry have 
also been discussed as future possibilities for Jefferson Lab, as has a more precise measurement 
of the M0ller asymmetry. The value of sin 2 6^(0) is particularly relevant for the interpretation 
of the parity violation experiments, since the vector coupling of the Z n Eq. (j3j) to electrons and 
protons is proportional to 1 — 4 sin 2 9w{.(fYf ~ 0.1 and is, therefore, highly sensitive to the value 
of the effective weak mixing angle. Indeed, as noted in Refs. [TUfTT] . one- loop contributions to 
sin 2 #ty(0)j ff reduce the magnitude of the M0ller asymmetry by roughly 40% from its tree- level 
value, an effect generated by the large logarithms discussed above. Moreover, the presence of 
these large logarithms in sin 2 #vy(0)j ff is universal to all low energy neutral current observables, 
though their net effect may be masked by other enhanced radiative corrections^]. Consequently, 
one would like to sum these universally enhanced contributions to all orders. Here, we do so 
using the renormalization group evolution for sin 2 #^(//) in the MS-scheme. 

The MS-definition of the weak mixing angle is, of course, not unique, and one may choose 
an alternate scheme in which to compute radiative corrections to electroweak observables. Nev- 
ertheless, the MS-scheme offers several advantages that motivate our adoption of it here. In 
particular, the MS-definition of sin 2 9w(^) follows closely the coupling-based definition in Eq. (JTJ) 
with a well-defined subtraction of singular terms arising in dimensional regularization, giving rise 
to expressions with a logarithmic //-dependence. This dependence is governed by a renormal- 
ization group equation (RGE), and choosing /i equal to the momentum transfer of the process 
under consideration will in general avoid spurious logarithms^]. As we discuss below, the evolu- 
tion of sin 2 ^(/ i ) can be related in a straightforward way to that of the QED coupling in 
the MS-scheme that has been thoroughly studied elsewhere. Doing so allows us to draw upon 
known results for the QED /3-function and - in conjunction with suitable matching conditions - 
to improve the precision of the Standard Model predictions for low-energy observables by incor- 
porating various higher order effects. Indeed, although the one- loop RGE for sin 2 9w{n>) has been 
well-studied by others, one emphasis of the present work is the inclusion of higher order QED 
and perturbative QCD contributions in a systematic way. We discuss the RGE in Sections [2] 

2 For example, the left-right asymmetry for parity-violating elastic ep scattering also contains large logarithms 
associated with the fermion anapole moment as well as non-logarithmic but large WW box graphs [12113114] . 

3 The presence of some logarithmically-enhanced radiative corrections, such as the anapole moment effects, 
cannot be cannot be eliminated by the RGE for sin 2 0vy (/■*)• 
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and [21 and matching conditions in Section HI 

A major complication arises when the contribution of the light quark flavors is considered for \x 
of order a hadronic scale, A, or smaller. In this regime, QCD corrections to the RGE (/3-function) 
cannot be obtained using perturbative methods. An analogous problem is well-known to arise 
for a(fj) when its value is desired at scales similar to or greater than A. We address this problem 
in Sections [5] and El as well as two appendices. We argue that the corresponding theoretical 
uncertainty is well below the anticipated experimental errors and provide a new estimate of 
this uncertainty that is substantially smaller in magnitude than the previously-quoted one [T5] . 
Non-logarithmic contributions and some more formal aspects are discussed in Section [7] while 
numerical results and a plot of sin 2 9 w (fi) for fi = \J\q 2 \ - along with conclusions - are presented 
in Section We note there our results may also be applied to the recent studies of deep- inelastic 
neutrino-nucleus scattering carried out by the NuTeV Collaboration [TB] . 

2 Leading order RGE analysis 

The quantity sin 2 9w(fy is related to s| = siia 2 9w(Mz) by, 

sin 2 § w (0) = k(0) sin 2 9 W {M Z ) = [1 + A/c(0)]s|, (7) 

where A/c(0) is a universal (flavor- independent) radiative correction. In this Section we are 
interested in logarithmic contributions to &(0) of the form, 

a . Ml 
In 



7r sin 2 9w m i ' 

and in particular in the scale that should be used in a and sin 2 9w, appropriate for re-summing 
the leading, large logarithms to all orders. These logarithms arise from scale dependent self- 
energy mixing diagrams where one external leg is a photon and the other one is a Z-boson. 
Thus, Vf acquires a compensating scale dependence, 

24^3/ HWhMriQi} in — - 2 



v f ( V/i 2 + A^)= $ f (jt) + -£LQ f YlW-lMtiQi] In - , 2 , (8) 



where = Nc = 3 (for quarks) and A| = 1 (for leptons) is the color factor. We have written 
the sum in general form to also allow chiral fermions and bosonic degrees of freedom (to put 
the W ± and Higgs ghosts on the same footing as the fermions and to facilitate the discussion in 
Section Ej), where the spin dependent factors ji are shown in Table [TJ With these conventions, 
the factor of 1/24 also appears in the lowest order QED /3-function coefficient, 



This implies the RGE, 



^0 = ^/E^7^Q„ (io) 
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field 


li 


real scalar 


1 


complex scalar 


2 


chiral fermion 


4 


Major ana fermion 


4 


Dirac fermion 


8 


massless gauge boson 


-22 



Table 1: Weight factors 7$ entering the leading RGE coefficients for the weak mixing angle. 



or in terms of the variable X = J2i NfjiViQi, 

dX a dp? 
X 247r y? 

where in the second equality we used Eq. (jUJ). This is solved by. 



da 
a 



sin 2 9 w (jj) 



sin 2 9 w (hq) + 



1 - 



d(/i ) 



&M vr ^ UJ ' EiW^Ql 

or using the explicit solution to the one-loop RGE in Eq. (Q we obtain the simpler form, 



sin 8 w (p) = sin 8 w (p c 



1 + 



247r sin 2 Owivo, 



53 N^iQi[Ti - Qi sin 2 6 w (n )} In 



f4 



;i2) 



(13) 



The result in Eq. (1130 re-sums all logarithms of 0(a n ln n ^) provided there is no particle 
threshold between \i and /jlq. To avoid reintroduction of spurious logarithms this solution must 
be applied successively from one particle threshold to the next. Crossing a threshold from above, 
the corresponding particle is integrated out, and one continues with an effective field theory 
without this particle. In contrast, changing from one effective field theory to another far away 
from the physical mass of the particle would not formally affect the truncated one-loop result, 
but it would spoil its re-summation. 

Eq. ([TBI applied to Ak(0) can be brought into the well-known form [1] 



A«(0) 



a 

7TS| 



1 /If 2 

gEW(T / -2g / 4)in^ 



/ 



/ 
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If. 
4 ' 



, M l 
In z 



Ml 



(14) 



where the sum is over all SM Dirac fermions excluding the top quark. The last term is the W ± 
contribution with its coefficient obtainable from Table [1] when a pair of massless gauge bosons 
(Tj = ±1) is combined with a complex scalar Goldstone degree of freedom (Tj = ±1/2). 



3 Higher order RGE analysis 

In this Section we will generalize the leading order analysis and re-sum next-to-leading (NL) log- 
arithms of 0(a n+1 ln n «a) and £>(«< ln n ^), as well as the NNL logarithms of 0(aa n s +1 ln n ^), 
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and the NNNL logarithms of 0(aa™ +2 In™ The leading order RGE (TITi) supplemented by 
terms of 0(a 2 ), 0(aa s ), 0(aa 2 ), and O(aa^) reads, 



dfi 2 



a 
24^ 




where in the case of quarks (n g is the effective number of quarks), 

; ' 2 '125 11 



7T" 



1 

55 



4^ 2 - + - + rf 



7T 



7T 



7T- 



10487 , , . 
1728" + 18 C(3) " 



/707 55 > , , . 
864 + 54« 3 >)"« 



n o + 

72 ' 
77 2 

— 72 

3888 q 



(15) 



(16) 



contains QED and QCD corrections [T7Tfl9] to the lowest order (non-singlet) vacuum polarization 
diagrams. For leptons only the term involving a is kept, while for bosons we restrict ourselves 
to the lowest order /3-functioii@, i.e. K w ± = 1. The second sum in Eq. (|15p is over Dirac quark 
fields, and 

a 3 

s 

71 3 



a 



216 9 v ; 



+ 0(a 



sJi 



(17) 



parametrizes the QCD singlet contribution. In a singlet (QCD annihilation) diagram two inde- 
pendent fermion loops are attached to the 7 and Z and connected to each other by gluons or 
photons. Due to Furry's theorem, connections containing a photon first arise at 0{a 2 a 2 s ) and 
can safely be neglected. Defining s 2 = sin 2 9w(n) we rewrite Eq. (fl5l) . 

/ \ 2 



, ds 2 
dfi 2 



a 

71 



KMQp 2 - TiQi) + as 2 (j2 Q Q 




Similarly, the RGE for at including higher orders reads, 



da 
dfi 2 



a 

7i 



24 



and we obtain, 



_d_ 

dfi 2 



24vr 



KiliTiQi 




(18) 



(19) 



(20) 



To facilitate the integration and to relate hadronic contributions as far as possible to the ones in 
a, we use Eq. ffl9|) again and eliminate all a s dependent terms in Eq. ffT6l) . With the coefficient^], 



Ai 



A 2 



\ E i^Ql - T iQi) = \T, n^x.q 2 - rm, (21) 



4 We do so because full two-loop electroweak calculations are generally incomplete and therefore only the leading 
order electroweak terms included in most current definitions of MS quantities. Moreover, the structure of the 
RGE would change relative to Eq. (|23p below, spoiling the corresponding solution ([25]) . Because the logarithms, 
]nMz/Mw, are not large, neglecting them in electroweak two-loop terms is numerically insignificant. 

5 The explicit factor of 1/2 in Ai compared to the coefficient in the leading order solution (fT2"|) arises because 
the electric charges in the denominator of the latter are summed over left and right chiralities while only left 
chiralities appear in the numerator. 
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energy range 


Ai 


A 2 


A 3 


A 4 


m < (j, 


9/20 


289/80 


14/55 


9/20 


Mw < fi < rh t 


21/44 


625/176 


6/11 


3/22 


rh b < fi < M w 


21/44 


15/22 


51/440 


3/22 


m T < fi < rhb 


9/20 


3/5 


2/19 


1/5 


rh c < fi < m T 


9/20 


2/5 


7/80 


1/5 


rh s < fi < rh c 


1/2 


1/2 


5/36 





rn d < /i < fn s 


9/20 


2/5 


13/110 


1/20 


m u < ji<m d 


3/8 


1/4 


3/40 





m M < fi < m u 


1/4 











m e < fi< rrif, 


1/4 











/i < m e 















Table 2: Coefficients entering the higher order RGE for the weak mixing angle. For the definition 
of quark threshold masses, rh q , see Section Below hadronic scales Ai is not defined through 
Eq. (12TT) and can be chosen arbitrarily; we chose the value, Ai = 1/4, to obtain A 2 = A3 = 0. 



and, 



A, 



ZiNfr&iQi-TiQi 



A, 



shown in Table [2] this can be brought into the form 



Ai £Q, 



_d_ 

dfi 2 



a 



X\ 3A3 



47T 



a 



A, 



E T * E^ 



hid + - = — . 

7T / 37T 



(22) 



(23) 



For the last term on the left hand side we have used the lowest order QCD /3-function coefficient 
and have defined 

iW .^i!hS +0(a3 . (24 ) 



7T 



2 36 33 - 2n n 



The solution to Eq. ( J23l is given by, 



a(iio) 



+ 



a(jj) 



7T 



A 2 , H 2 , 3A 3 a(fj) , , 
3 //g 4 a{fio) 



■ (25) 



Hadronic uncertainties are induced through d(/x), through the relative values of the light quark 
threshold masses, m u , m,j, and m s (these are needed as they determine the change of the coeffi- 
cients Aj according to TableEJ), and through the singlet contribution proportional to a([io) — a(fi). 



4 Matching conditions 

At the threshold of fermion / we find, 

• 2n I \- Ol(lTli) 

sin Vw(mi) 



&(mi)- 



2Q1 



1 - 



a(mi)- 



(26) 



where the plus (minus) superscript denote the effective theory including (excluding) fermion /. 
While the one and two-loop /3-function coefficients are well-known to be renormalization scheme 
independent, the matching conditions are renormalization scheme and even regularization scheme 
dependent. The MS-scheme is defined by dimensional regularization which generates no 0(a) 
matching terms for scalars; with the usual additional requirement^ that the Clifford algebra is 
kept in four dimensions the same holds for spin-1/2 fermions. We include the RGE matching 
conditions for a at the threshold for fermion / at the orders a 2 , aa s 



1 



Q) 



7T 



15 
16' 



7T 



(27) 



(Nf-l)&.(m f ) 



7T 



13 & s (m f ) /655 _ 3847 
12 7i [ 144^ ' 864 



361 295 E q ^f Q\ 

-n a + 



1296' q ' 1296 Q} 

Here, n q is the number of quarks including the threshold quark© /. Eq. (1261) will then induce the 
corresponding matching contributions to the weak mixing angle. 

In contrast to fermions and scalars, gauge bosons induce an 0(a) threshold shift 



a. 



C{R) 
12tt : 



f28l 



where C(R) is the quadratic Casimir of the (in general reducible) gauge boson representation, 
normalized such that, e.g., C( adjoint [577 (iV)]) = N. Thus, integrating out the W ± bosons 
induces a shift© in the electromagnetic coupling, 



1 

67T 



a(M w )+ a(M w ) 
and (generalizing Eq. (1261) appropriately) in the weak mixing angle, 

. a{M w ) + 



(29) 



sin 2 6 W (M ] 



a(M w )- 



cos 2 9 W (M W ) 



(30) 



5 Hadronic contribution 

The ambiguity in the values of the three light quark masses plus the singlet contribution to 
Eq. (|25|) introduce four sources of hadronic uncertainties in sin 2 9^(0). This problem is familiar 

6 Othcr definitions do occur in the literature, however. 

7 In writing this equation we assume that mj is an MS-mass (which is free of renormalon ambiguities and 
assures a better convergence of the perturbative series) to the extent to which QCD effects are concerned, but 
a pole mass for both leptons and quarks with respect to QED (to comply with standard conventions in the 
literature). This results in a somewhat awkward definition for quarks but is of no importance in practice since 
the 0(a 2 ) corrections are very small. 

8 This shift is an artifact of using modified minimal subtraction in dimensional regularization. It is precisely 
canceled against a conversion constant 23 24J which relates the MS-scheme to the DR-scheme. The latter is 
defined by dimensional reduction and is used in supersymmetric theories. 
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from the evaluation of a(Mz), where it can be addressed by relating it via a dispersion relation to 
e + e~ — > hadrons cross section data, or (using in addition isospin symmetry) to r decay spectral 
functions. The same strategy could in principle be applied here, except that the experimental 
information would have to be separated in charge 2/3 (u) vs. charge —1/3 (d and s) quarks, or 
assuming isospin symmetry, in s vs. the first generation quarks. This is a difficult task: e.g., 
a K ± final state in e + e~ annihilation can be produced directly through an s7 M s current or by 
splitting of a gluon radiated off a quark originating from a vry^u + dj^d (isoscalar) current. It 
would be even more difficult to isolate the singlet contribution. As for r~ decays, the relatively 
large m s induces a sizable axial- vector contribution to final states with strangeness, S = — 1. At 
least presently, it cannot be cleanly separated from the vector contribution which is the relevant 
one for the problem at hand. We will therefore follow a different strategy and show that the four 
uncertainties can be traded for the uncertainties associated with (i) the value of a(Mz), which 
we denote S a sin 2 8w(0), with (ii) the separation of strange and first generation quark effects, 
indicated by 5 S sin 2 6^(0), with (iii) deviations from isospin symmetry, <5cvc sin 2 6W(0), and with 
(iv) Zweig (OZI) rule deviations, 5ozi sin 2 #w(0). 

As discussed in Section [3j the contributions of leptons and heavy quarks can be computed 
unambiguously in perturbation theory. Indeed, using QCD MS massesj, rh q (m q ), will provide a 
small truncation error, #pqcd sin 2 §w(0), in the perturbative expansion (cf. the small matching 
coefficients in Eq. ( 1271) ). Moreover, the numerical values of rh q (rh q ) can be determined to sufficient 
precision that is not limited by uncertainties of the order of hadronic scales. Therefore, we 
compute first sin 2 0w(p) from sin 2 9w(Mz) (which can be taken from Z pole experiments), where 
p, corresponds to a scale where the heavy flavors (b, c, and r) are integrated out, and at which 
we still have sufficient confidence in the convergence of perturbative QCD (i.e., of order 1 GeV). 

We constrain the contributions of light quarks to sin 2 dwifi) phenomenologically. Our strategy 
is to employ the u, d, and s quark contributions to ce(fi), Ace^(p,), as a constraint, and to find 
upper and lower bounds on the strange quark contribution relative to the contributions of the 
first generation quarks. In the following we assume isospin symmetry and a vanishing singlet 
contribution. Deviations from these assumptions will be addressed in Section [6j 

To facilitate the discussion, we will adopt definitions of threshold quark masses, m u , rh d , 
and m s , such that Eq. (|25|) remains valid with trivial matching conditions, af(rh q ) = a~(fh q ). 
Thus, fhu, rhd, an d fh s define the ranges in Table [2j Their values can be constrained phenomeno- 
logicalljo, but their relation to other mass definitions, such as constituent masses or current 
masses, cannot be written down in a perturbative sense. One combination of the three light 
quark threshold masses is constrained to reproduce a(jl)/a>. If we assume isospin symmetry, 
m„ = rhd, and a vanishing singlet contribution, we have only one unknown parameter, say rh s , 
to describe S s sin 2 6^(0). Before we use physics arguments to constrain fh s , we compute rh c and 
fhb perturbatively to gauge the behavior of heavy quarks. To order a 2 we have [20f2~T] . 



In 



| & s (fx) 
m 2 (/i) 7r 



13 ji 



2 



In 



12 m 2 (/i) 



a 

+ 



2 



n 2 



655 3847 361 295 E^/Q 



rC(3)-^T- + n g — - + 



144 s v ' 864 q 1296 1296 Q 2 f 



f 



0, (31) 



9 See Ref. 1251 for an earlier determination. 
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which implies, 



m = m{m) exp 



13 &.(m) &l 3847 361 295 E^/ Q q 

' 24 a s (m) + n 288^ [ bbb ^ 6 + 9 9 + 9 



(32) 



With the input values (obtained from a global fit to precision data), a s [Mz) — 0.1214 ± 0.0018, 
m c (m c ) = 1.285±g;g|?, and m b (m b ) = 4.205±0.031, we find m c = 1.176 GeV and m b = 3.995 GeV. 

The heavy rh s limit 

To obtain a lower limit on the strange quark contribution to a(/i), we consider the case in which 
the strange quark is assumed to behave like a heavy quark. In this case, fh s would be related to 
Mfj, in a similar way as rh c (or m c (m c )) is related^ to Mj/^. Defining ^ q = 2rh q /M ls where M 1S 
is the mass of the IS qq resonance, we have that asymptotically £ q — > 1 for fn q — > oo and £ g — > 
for fh q — > (in the chiral limit the quark contribution is logarithmically divergent). Thus, for a 
heavy quark ^ q ~ 1, while for a light quark ^ < 1. Also, we expect fh\ < rfi2 ==>• £i < £2- As 
an illustration, with the numerical values of rh c and rh b from above we obtain £ c = 0.759, £& = 
0.845, and m s = £ s M^,/2 < ^M^/2 = 387 MeV. As a refinement we introduce scale dependent 
QCD correction factors, Kq CB = KQ CD (p,,fh q ), where -Kqcd(/^i, Jh) denotes the average QCD 
correction to the QED /3-functions for RGE running between scales fx\ and Jii- One thus expects 
rhi < 1^2 =>• -Kqcd > -Kqcd- Since Eq. (132]) applied to m c still shows satisfactory convergence, 
we can safely choose p, = m c , 



A s a(m c ) = Qj£tf&D ^ g > Q^qcd In g > f^ CD In ^ = 6.9 x 10" 4 . (33) 



s 



For the numerical evaluation we have used the QCD correction in Eq. (fT6l) applied to the effective 
theory with n q = 3 quarks, 

K s ^ K c -I i «jW , 103 &g(m c ) 1979 gg(mj 

^qcd > QCD = 1 + + ^ + "576 ^3 = i' 209 ' (34) 

and we have used a s {fh c ) = 0.469, again corresponding to a s (M z ) = 0.1214. 

The SU(3) limit 

Since m s > nid > m u at any scale and in any reasonable definition^ and scheme, we conclude 
that the SU(3) symmetric case, rh u = rhd = rh s , implies an upper limit on the relative strange 
quark contribution to a(£L), 

A s a{p) < \ [A u a(fi) + A d a{fl) + A s a{fl)) = ^Ad (3) (/i). (35) 
6 6 



10 More precisely, QCD sum rules relate fh c (nib) rigorously to a weighted sum over VE" (T) resonances plus a 
continuum contribution. For the present consideration we restrict ourselves to the lowest lying resonance which 
carries the largest weight. 

11 This statement holds because small non- universal mass renormalization corrections from QED and the elec- 
troweak interactions can be neglected. 
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This crude limit can be strengthened by considering the phenomenological constraint, 



2a 

71 



(Ql + Q't)K$ D \n 



2p 



+ Q 2 S K S QCD In 



2p 



Ad( 3 )(/2), 



(36) 



and by imposing SU(3) symmetry through £ u — £ d — £ s and Kq CD = Kq CD = Kq CD . This 
maximizes the ratio of the strange quark contribution to the one of the first generation quarks, 



Ql 



A u+d a(jl) ~ Ql + Q 2 d 



1 + 



QCD m c -^ QCD 111 M 



In M± — (fC u > d _ K 



s \ In 2jJ 



QCD m t a jM- tt 



(37) 



All three corrections terms to the charge square ratio are indeed negative, and we used the second 
one as an improvement. In the limit £ s — > 0, Eq. fl36l) would reproduce relation fl35|) but now we 
have the constraint, 



2/2 



M* /6 Ml /6 



exp 



3?rAa( 3 )(/i) 



4 aK 



QCD 



(3* 



For the numerical evaluation we convert the contribution to the on-shell definition of a(Mz), 
Aa haxi (M z ) = 0.00577±0.00010 [271128] . from the energy range up to 1.8 GeV, to the MS-scheme, 



Ad< 3 >(1.8 GeV) = Aa had (M z ) + ^ 



a 



55 



o+ - + T- 77T- 4 C(3) + 



2m? 



7T 



4vr / V 12 



1.8 GeV) 2 



+ 



a] ( 34525 9 A/X 715 25 A/X „, A A , 

C(2) C(3)H C(5) +F(m c ,m b 

864 4 SV ; 18 sv ; 3 sv ; v ' ' 



7T 



I" 



0.00831 ±0.00010, 



where all MS running couplings and masses are to be taken at // = 1.8 GeV. F(rh c ,rhb) contains 
decoupling charm and bottom mass effects [T9|2 9 30 j. We choose again p, = m c and use the 
4-loop RGE to obtain, 



Ad (3) i 



0.00678 ±0.00010, 



(39) 



which using the SU(3) bound (1381) corresponds to £ s > 0.470 and m s > 240 MeV. Inserting 
Eq. (J3HD into Eq. ([36]) in the SU(3) limit yield@ 



A.a(m c ) < 



Aq( 3 >i 



5 a 



6 



27 7T QCD M (JJ 6 27 tt QCD M 



9.9 x 10" 



(40) 



We have assumed ideal u — (j) mixing, i.e. that the <fi resonance is a pure ss state. Allowing a 
non-ideal mixing angle, e = 0.0548 ± 0.0024 ^ (see Appendix IA1 . shifts the masses and 
to be used in Eqs. fl3"3"|) and (1401 by less than 1 MeV and yields a negligible effect. 



12 Exact SU(3) symmetry would imply M p — M u ; since we arc interested in an upper limit on the strange quark 
contribution we have choose the (larger) phenomenological value of M w . 
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Implications 

From Eqs. (133]) and (HOI) we conclude for the strange quark, 

A s a(m c ) = (8.4 ± 1.5) x 10~ 4 , 
and for the light quarks, 

Aa {2 \m s ) = A« (3) (m c ) - 6A s d(m c ) = 0.00172 =f 0.00090, 



m s = 305+^ MeV. (41) 



m d = rh u = 176 ± 9 MeV. (42) 



These results can be used in the master equation (1251) . As an illustration, the SU(3) symmetric 
piece is well approximated by, 



s (m s ) - s (m c ) 



s 2 (m c ) 



6A s a(m c ) 



2a(m s 
3vr 



1 



s 2 (m r 



h^ l + W, (43) 
m 2 . \ 47r J 



which is obtained with the help of Table [21 and where the last term is the leptonic (e and //) 
contribution. Similarly, the 577(3) breaking piece reads, 



s 2 (m d )-s 2 (m s ) 



9 

20 



s 2 (m x 



Aa^i 



2a(m d ) 



3tt 



s 2 (m. 



47T 



(44) 



where we neglected the singlet piece involving a. 



6 Uncertainties 



From Eq. (159]) . as well as the sum of Eqs. (|4"3"]1 and (144]) . we can bound the uncertainty induced 
by Aa' 3 '(m c ), 

■1 



5„ sin 



--s [m c ] 



±3 x 10 



-5 



(45) 



Similarly, from the first Eq. ( |4"2]) and from the comparison of the coefficients in Eqs. ( 1431) and 
(l4"4"j) . we can estimate the uncertainty induced by Aa' 2 '(m s ), 



<5 s sin 2 M0) w ^-5Aa (2) (m s ) = ±5 x 10" 5 . 



(46) 



The singlet contribution 

We obtained the theoretical bounds on rh s that are the basis for the error estimate (T46]) by assum- 
ing isopsin symmetry and a vanishing singlet contribution. We now relax the latter assumption 
and allow OZI rule [3TH33] violation which leads to processes such as <p 7r °7 decays and which 
translates on a diagrammatic level to QCD annihilation (singlet) topologies. For charm and third 
generation quarks singlet contributions are tiny and easily included using Eq. (ITT]) or Eq. (124]) . 
We can then proceed with the effective theory containing only the three light quarks. Notice, 
that due to Q u + Q d + Q s = 0, there is no singlet contribution (to the /^-functions of neither a 
nor sin 2 9^) in the limit of exact SU(3) symmetry (see also the sixth entry for A4 in Table [2]). 
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Moreover, allowing SU(3) breaking effects — but still working in the isospin symmetric limit - 
will not directly affect the RGE (120"]) because T u + T d = 0. The explicit singlet piece in Eq. ([25]) 
is an artifact of employing Eq. (TT9"]) and cancels the implicit singlet piece contained in the term 
proportional to Ai. Not being able to isolate the implicit piece phenomenologically or to calculate 
the explicit piece in the non-perturbative domain introduces an additional uncertainty, which we 
now argue is rather small. 

Perturbation theory provides an order of magnitude estimate if one assumes that the leading 
order perturbative coefficient is of typical size and not accidentally small. Then one would find 
for the singlet contribution, 



<$ozi sin 2 6^(0) ~ A 4 - 

7T 



3^ 

7T 



2 _5_ 11^2403) 
324 33 - 2n g ' v 1 



where the QCD expansion parameter in square brackets has been assumed to have grown in the 
non-perturbative regime to a number of 0(1), and where n q = 2 and A4 = 1/20 correspond to 
the effective field theory with the strange quark integrated out. More generally, based on results 
of Ref. [3l] we anticipate that in leading order in Nc the singlet terms are of the form, 



<$ozisin 2 6W(0) ~ A 4 - 

7T 



7T 



T F C F N C „ lot 5 



where the QCD group factors are T F = 1/2, Cp = 4/3, Ca = N c = 3, and where the coefficients 
C n are expected to be of 0(1). An alternative form can be written down relative to Aa^(m s ) 
in Eq. (02]), 

5ozi sin 2 M0) ~ ^q^qi ^T = 4 Aa(2)( ^ s) = (2 ' 6 T lA) X 10_5 ' (49) 

which incidentally gives the same result. These forms exhibits all QED charges and leading QCD 
group factors explicitly, which combined lead to a suppression of the singlet contribution by two 
orders of magnitude relative to the non-singlet contribution. Thus, the smallness of the estimate 
in Eq. (|47|) is in part due to Ci ~ —0.043 1 (which may or may not reflect the typical size of 
the other C n ), and in part due to the suppression factors displayed in Eqs. (14"51) and (|4"5|) which 
will apply at any order. In particular, OZI rule violating effects are absent in leading order in 
the 1/Nc expansion. 

In Appendix |B] we will test these order of magnitude estimates by studying the masses and 
mixings of vector mesons (which strongly dominate the real parts of the vector current correla- 
tors). The results obtained there turn out to be in line with the estimate (T4T1) . but conservatively 
we base our final uncertainty on Eqs. (1481) and (1491) and take, 



5ozi sin 2 9 W (0) = ±3 x 10" 5 . (50) 



Isospin breaking 



So far we have assumed exact isopsin symmetry. Recall that the SU(3) limit serves to maximize 
the RGE running of sin 2 6w by minimizing the effective down-type quark masses relative to rh u . 
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Allowing fhd > m u can therefore only strengthen this limit. Thus, the uncertainty associated 
with isospin symmetry breaking, £cvc sin 2 9w(0), is asymmetric. 

As for the heavy fh s limit, we proceed by considering the hypothetical reference case, fhd = 
ffis 7^ Tfiu- In our framework this corresponds to maximal SU(2) (CVC) violation, i.e., SU(2) 
breaking is of the same size as SU(3) breaking. The inequality, A.a <y2 \m s ) < 0.00262 (see 
Eq. ( BSD ), would be replaced by, 

Aa w (m d ) < 0.00262, (51) 

which bounds the up quark contribution for energy scales below the down-type quark effective 
masses. Ai = 3/8 now to be used in Eq. (T25]) in place of Ai = 9/20 in the isospin symmetric case 
would cause a shift, 

<5 CV csin 2 M0) = -^u W (™ d ) > -2 x 10~ 4 . (52) 
A measure of SU(2) breaking relative to SU(3) breaking is given by the ratio, 

M 2 K , ± - M 2 K , 



M\, ± - M 2 
which leads to the estimate, 



0.04, (53) 



5 CV csin 2 M0)=^ xi(T 6 , (54) 
and shows that isospin breaking affects our analysis at a very small level. 



Other uncertainties 

In the perturbative regime we used theory in place of experimental data, which induces two 
kinds of uncertainties: purely theoretical ones and parametric ones from the input quark masses 
and a s . The former include the errors associated with the truncation of perturbation theory 
and with non-perturbative effects. We estimate their size to about ±7 x 10~ 5 . However, this 
uncertainty is already included in Eq. (1391) where it propagates properly correlated to the error 
of the low energy weak mixing angle. The uncertainties in the quark masses induce an error 
of ±4 x 10~ 5 which is dominated by m c (m c ). The uncertainty in a s induces an error of the 
same size. In practice, these parametric uncertainties and the one from Eq. (1391) are included in 
global fits to all data where these parameters are allowed to float subject to experimental and 
theoretical constraints, and where correlations are naturally accounted for. The same applies to 
the experimental uncertainties in sin 2 6y/(Mz) and Mz- These have almost no effect on k(0), 
but if the absolute normalization of s 2 (0) is required they induce errors of ±1.4 x 10~ 4 and 
±1.4 x 10~ 5 , respectively. 

The theoretical uncertainties (T45|) . (T46|) . (|50l) . and (|54j) added in quadrature yield a total 
theory error, 

<W y sin 2 M0) = ±7x 10- 5 . (55) 

This is almost an order of magnitude more precise than the result obtained some time ago in 
Ref. [15j. Using our results in a global fit to precision data yields, 

sin 2 § w (0) = 0.23867 ± 0.00016. (56) 
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The central value coincides with Ref. [36] where a seemingly independent definition of the low 
energy mixing angle (based on gauge invariance and the pinch-technique) is introduced. We will 
comment more on the relation between our work and Ref. [36J in the following Section. The 
uncertainty in Eq. (156]) is completely dominated by the experimental uncertainty, 

6 sin 2 9 W (M Z ) = ±1.5 x 1(T 4 . (57) 



7 Other considerations 



Some time ago, the authors of Refs. [TOlfTT] suggested that the use of an appropriate, scale- 
dependent effective weak mixing angle could provide a useful means of comparing the results 
of various neutral current experiments at the Z°-pole and below. By now, it is conventional to 
compare the value of an effective weak mixing angle extracted from experimental results with its 
predicted value in the Standard Model (see, e.g., Ref. [8]). More recently, it was observed [36] 
that the effective weak mixing angle derived from the sum of Z-7 mixing diagrams evaluated in 
the gauge and frequently used to interpret low energy neutral current experiments is not gauge 
invariant. This sin 2 9(q 2 ) eS is defined analogously to Eq. ([5]) with a scheme- and g 2 -dependent 
form factor, k(q 2 ), 



sm 



2\ eff 



Re[% 2 , n = M z )\ sin 2 § w (fi = M z ) 



(58) 



In the MS-scheme, both k and sin 2 9w depend on the renormalization scale fi, while k also carries 
a (^-dependence. The authors of Ref. [36j note that the k(q 2 ,fi) form factor naively-defined in 
terms of Z-7 mixing depends on the choice of electroweak gauge so that the corresponding 
sin 2 9(q 2 ) eff is not a physically meaningful quantity. By itself, this gauge dependence is not 
particularly problematic, since for any physical observable — such as the parity violating M0ller 
asymmetry computed in Ref. [10] - it is canceled by the gauge dependence of other radiative 
corrections, leaving a gauge independent result. Nevertheless, if one wishes to isolate a particular 
class of radiative corrections, such as those entering sin 2 9(q 2 ) eff , one usually prefers to discuss 
gauge invariant quantities, especially when the comparison of different experimental results is 
involved. 

The authors of Ref. [36] show that one may, indeed, obtain a gauge independent k(q 2 ,n) 
by including the so-called "pinch parts" of various one-loop vertex and box diagrams that are 
process independent and that compensate for the gauge dependence of the naive k form factor. 
Here, we comment on the relationship between sin 2 9(q 2 ) eff of Ref. [36] and sin 2 9w(n) discussed 
in our work and observe that they are identical at one-loop order. 

For \q 2 \ < Myy, the gauge invariant form factor k of Ref. [36J is given by, 



k(q 2 ,fi) 



PT 



1 + 



a 



2ixs 2 z 



, M 2 

n 2 



7 



-,Y. N fQA T f- 2 QfS 2 z) + oC Z + - 



a 



2vrs| 



2j2N c f Q f (T f -2Q f s 2 z )I f (q 2 ) 
f 



7.2 Mi M z 



r 2 
+ — 

3 



(59) 
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where, 

m 2 — q 2 x(l — x) 



I f (q 2 ) = ^ dxx{\ - x) In 7 M , (60) 

and where the PT superscript in Eq. (|59[) indicates the gauge invariant "pinch-technique" defini- 
tion of the form factor. For \i = M z , the second term on the right hand side of Eq. (1591) vanishes. 
The integrals If{q 2 ) in the third term generate the large logarithms containing fermion masses 
that one would like to re-sum. As discussed in Section [TJ this re-summation is accomplished by 
choosing \i ~ \J\q 2 \ rather than /i = M z , thereby eliminating these logarithms from k(q 2 ,fi) PT 

altogether and moving them instead into sin 2 9w{^)i which we analyze in this paper. The RGE 
for sin 2 #jy(/-0 then provides for the desired re-summation. Similarly, the term proportional to 
In M\jM w in Eq. (I59p corresponds to the weak gauge sector contributions to the RGE running 
from fj, = M z down to M w . Below this scale, the heavy gauge bosons are to be integrated out. 

It is not too surprising that the logarithms appearing, for example, in Eq. (fT4"j) . are identical to 
those obtained from the PT since it has been shown [37] that the asymptotic behavior of effective 
coupling constants directly constructed from the PT self-energies are automatically governed by 
the renormalization group. Now we observe that even the non-logarithmic piece in the third term 
of Eq. (|59|) can be understood in the context of the renormalization group, except that in this 
case it arises from RGE matching rather than RGE running. In Ref. [36] the c|/3 term results 
from combining the pinch parts of the one-loop vertex and box graphs with the remaining weak 
gauge- dependent contributions to the 7-Z'-mixing tensor. The precise value for this /i-independent 
constant follows from the requirement that k(q 2 ,fi) PT be gauge invariant. In our treatment of 
the running sin 2 9w(n), this same constant is generated by the threshold corrections at \i = Mw 
given in Eqs. (1291 and (l30l . Indeed, use of the RGE with appropriate matching conditions may 
provide a more direct route for obtaining the results of Ref. [36] while allowing one to generalize 
it to include various higher-order effects as we have done. 

It follows as a corollary that the PT applied within the DR-scheme (compare the last footnote 
in Section H]) should not yield any constant terms at one- loop order. As a particular application, 
one may consider the correspondence between the two treatments at // = 0. Eqs. (fT2"l) . (1291) . 
and (1301) show that the relation between sin 

6 W {0) and sin 2 9 W {M Z ) 
sm 2 § w (0) = sm 2 § w (M z ) (61) 



a 

+ - 

7T 



- V N c f QAT f - 2Q f s 2 z ) In - (- - -s 2 z ) In ^ - -c| 
6 y /VA S V/ z; m) V24 4 z ) M w 6 z 



is the electroweak analog of the relation between the electromagnetic fine structure constant, 
a = a(0) and a(M z ). Note that a similar relation would hold for other definitions of the weak 
mixing angle and the corresponding definition of the running QED coupling in the same scheme. 
For example, different conventions for the treatment of heavy top quark effects [5] would affect the 
definitions of sin 2 9w(M z ) and a(M z ), but the right-hand side of relation (|6T|) would also have 
to be modified with the net effect that the left-hand side would remain unchanged. Nonetheless, 
we reiterate that the definition (|6T|) is gauge invariant because it agrees with sin 2 #(g 2 = 0) eff , 
and that the analysis of Section [3] has allowed us to incorporate higher order effects in a and 
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a s into sin 2 9 W . Note also that if sin 2 9w([i) or sin 2 9(q 2 ) eff are used in low \q 2 \ amplitudes, care 
must be taken to consistently include other radiative corrections in the same scheme. 

While our study has focused on sin 2 9w{n) appropriate for low \q 2 \ processes, it is also worth 
commenting on the running of the weak mixing angle above the weak scale. For \q 2 \ ^> M§ , it is 
most appropriate to work in a basis involving the "primordial" SU(2)l and U(l)y gauge bosons 
and the /3-functions for g and g' . Starting from Eq. (CQ) one obtains the RGE, 

§2 da _ & d^ = b 1& 

dt dt IT 4^ 7T 2 ■? ' v ; 

J 

where t — In //, and where 62 and &2j are, respectively, the one- and two-loop /3-function coefficients 
involving SU{2) L (see, e.g., Ref. [38]). The solution to Eq. (|62|) can be written in the same form 
as Eq. (125]) . We note, however, that a naive application of the RGE ( 162]) to scales /1 -C M w 
would not re-sum all the large logarithms associated with the low \q 2 \ radiative corrections. For 
example, from Eq. ( 1621 one obtains, 

y t 2 

So long as both members of a quark doublet are included in the effective theory, this result is 
equivalent to the expression in Eq. (I2"T]) . since Q = T + Y and Tr (TY) = 0, where Y denotes 
hypercharge. However, for \i lying between the masses of two doublet members, this equivalence 
no longer holds, and only Eq. (T2~T1) will lead to a full re-summation of the large logarithms. 



8 Conclusions and outlook 

With the completion of the precision electroweak programs at LEP 1, SLC, and LEP 2, precision 
measurements of low energy neutral current observables have taken on added interest in recent 
years. A useful way to compare the results from existing and prospective experiments is to 
extract the value of the weak mixing angle that each would imply, assuming no other physics 
than that of the SM. The extent of their agreement with the SM prediction for this quantity 
provides important information about both the SM as well as the various scenarios that might 
extend it. 

The impact of this low energy precision program depends on both the precision of the var- 
ious experiments as well as that of the SM predictions. In this study, we have attempted to 
refine the latter by giving the appropriate low energy running weak mixing angle in the MS- 
scheme, sin 2 6*iy (/-*■)• By using this quantity and taking /i 2 ~ \q 2 \, one is able to re-sum various 
logarithmically-enhanced contributions that would otherwise appear in the radiative corrections 
for \q 2 \ Myy, thereby reducing the truncation error associated with the perturbative expan- 
sion. At one- loop order, this re-summation reproduces the result of Ref. [36] , but we have been 
able to generalize that work to include higher order contributions in a and a s . We have also 
provided an extensive analysis of the non-perturbative hadronic contributions to sin 2 9w(fi) for 
H ~ and argued that the associated uncertainties enter below the 10~ 4 level. 

The resulting scale-dependence of sin 2 9-wip) for fi = \J\q 2 \ with q 2 being the four-momentum 
transfer squared is shown in Fig. [H The various discontinuities in the curve correspond to the 
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Weak Mixing Angle 



Scale dependence in MS scheme including higher orders 




Q [GeV] 

Figure 1: Scale dependence of the weak mixing angle in the MS renormalization scheme. 

thresholds discussed above, while the size of the theoretical uncertainty in the curve corresponds 
to its thickness. 

As a particular application, we obtained a definition of the mixing angle in the Thomson limit, 
sin 2 9w = sin 2 #vk(0) whose relation with the value determined at the Z-pole sin 2 9w(Mz) is the 
electroweak analog of the relation between the fine structure a = &(0) constant and aiMz)- This 
definition also coincides with the gauge invariant definition recently constructed in Ref. [36J, and 
its numerical value is 

sin 2 9 W = 0.23867 ± 0.00016, (64) 

where the error is dominated by the experimental error from Z pole measurements. From its 
relation to the MS-scheme mixing angle at the Z scale, sin 2 9w(M z ), by Eq. ©, we obtain 

Ak(0) = 0.03232 ± 0.00029, (65) 
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where now the error is purely theoretical. Finally, using the relation [35] between sin 2 §w(Mz) 
and the effective leptonic mixing angle, sin 2 #y ff , defined in Eq. (JSJ), we obtain, 

sin 2 6 W - sin 2 6 e cS = 0.00718 ± 0.00007. (66) 

The error in this relation (which is the main result of this work) is an order of magnitude below 
current and anticipated experimental uncertainties and considerably smaller than the uncertainty 
quoted in Refs. [T5] . 

To illustrate the impact of this result on particular observables, we consider the weak charge 
of the proton, Qwip), that will be determined using parity- violating elastic ep scattering at 
the Jefferson Lab. Up-dating the recent analysis of Ref. [39], we obtain the Standard Model 
prediction 

Q w (p) = 0.0713 ± 0.0006 (input) ± 0.0003 (Ak) ± 0.0005 (Z7 box) ± 0.0001 (WW box) 

= 0.0713 ±0.0008, (67) 

where the experimental uncertainties in sin 2 §w(Mz) and m t ("input"), the theoretical hadronic 
uncertainties in Ak(0) and the Z'-f box graph, and the uncertainty from unknown higher-order 
perturbative QCD contributions to the WW box graphs are shown separately and combined 
in quadrature in the second line. Use of the previous estimate for the uncertainty in Ak(0) of 
±0.0025 would lead to a theoretical (total) error in the Qw(p) prediction roughly five (three) 
times larger than in Eq. ( 1671) . and in this case, one could neglect the uncertainties associated with 
other radiative corrections. In light of our analysis, however, the uncertainty in the Standard 
Model prediction is now three and a half times smaller than the anticipated experimental error, 
and theoretical uncertainties associated with hadronic contributions to other radiative corrections 
become relatively more important. 

In the same vein, the interpretation of the prospective parity violating deep inelastic measure- 
ments will require a careful analysis of higher twist and isospin breaking corrections, especially 
given that the latter may be responsible for the present discrepancy between the NuTeV result 
for sin 2 W (// ~ 3 GeV) and the SM prediction. Our analysis applies to the deep inelastic regime, 
as well, where due to the higher energies involved, the uncertainty in sin 2 6 l iy(/i) is even much 
smaller since no complications from non-perturbative contributions arise in this case. Given the 
level of experimental effort required to carry out these precise low energy measurements, per- 
forming a theoretical analysis of these effects at the level we have attempted to do here for the 
weak mixing angle seems well- worth the effort. 
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A Updated value for the mixing angle e 

Here we determine the phenomenological value of the u-<f> mixing angle, e, defined by, 

\uu) + \dd) \uu) + \dd) 

\u) = cose -= sin e | ss), \<p) = sine -= h cose|ss), (68) 

\/2 \/2 

where e = is referred to as ideal mixing. One way to extract it is to use £77(3) flavor symmetry 
and first order breaking applied to the vector meson octet mass spectrum. An advantage of 
this method is that it can be calibrated against the ground state baryon octet, for which Fermi 
statistics precludes the mixing with an SU(3) singlet state. The mass of the £77(2) singlet, M\, 
is therefore predicted in terms of the masses of the other electrically neutral octet members, 

M A = -(2M„ + 2M=o - M s o) = 1105.4 MeV, (69) 
3 



which reproduces the experimental value |40j, Ma = 1115.7 MeV, within 0.93%. Analogously, 
the Gell-Mann-Okubo mass formula [4Tf4*2] yields the octet-octet component of the mass matrix 
for the isosinglet ground state vector mesons, 

M 8 2 8 = -(4M^ - M 2 ) = (933.69 MeV) 2 X [1 ± 0.0008 ± 0.0020 ± 0.0121 ± 0.0093]. (70) 

The first error is from the experimental uncertainty in the masses, which are taken from Ref. [40] 
except for the mass, M p o = 775.74 ± 0.65 MeV, and the width, T p o = 145.3 ± 1.4 MeV, of the p° 
resonance for which we averaged Refs. [4314*4] . The broadness of some of the resonances involved 
introduces an ambiguity as for what definition of mass one should use in the Gell-Mann-Okubo 
formula. For the central value we have chosen the peak position, 

where M and T correspond to the usual definition of a relativistic Breit-Wigner resonance form 
with an s-dependent width, i.e. the one used by the Particle Data Group [40]. The second error 
in Eq. f[70l reflects the size of the shift obtained by replacing M by M. The third error is due to 
possible isospin breaking effects which we estimated by using M\,± in place of M|-, . The last 
error quantifies the limitation of the method and is given by the calibration against the baryon 
octet as discussed above. Adding all errors in quadrature, we obtain for the SU(3) octet-singlet 
mixing angle [40], 9y, 

tan 2 ^ = ^^= 0.646;^, (72) 
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which translates into the two solutions 13 !. e\ = 0.061 =F 0.032 and e 2 = —1.292 ± 0.032. Alterna- 
tively, one can compare the branching ratios of the u and <fi resonances decaying into 7r°7 



tan 6 ~ M 3 \Ml - M* ) B{u - 7T° T ) T w " (3 - 01 ± °- 26) X 10 ' (73) 
which gives, |e| = 0.0548±0.0024. Comparison with the previous method singles out the solution, 

e = +0.0548 ± 0.0024. (74) 



The sign and magnitude are also consistent with various other methods (46] and the previous 
analysis in Ref. [45] . 



B Phenomenological approach to OZI suppression 

In Section [6] we argued that the OZI rule can at least partly be understood as a result of 
group theoretical suppression factors relative to OZI allowed processes. Using the result of 
Appendix [A], we now wish to study OZI rule violating contributions to the mass matrix of 
ground state vector mesons. These mesons dominate the electromagnetic current correlator 
at hadronic scales, and should therefore serve as a means to quantify OZI rule suppressions 
phenomeno logically. Throughout this Appendix we work in the isospin symmetric limit. 

In the flavor basis, (\uu), \dd), \ss)), we write the mass matrix in a form which is similar to 
the one discussed in Ref. [47], 

M= I B A + B B + C I. (75) 




The parameters A and B respect SU(3) symmetry, while C and D break it. The off-diagonal ele- 
ments, B and C, parametrize flavor transitions, and can only be generated by QCD annihilation 
(singlet) diagrams. The parameter D receives contributions from the strange quark mass, as well 
as dynamical contributions of both singlet and non-singlet type. The difference to Ref. [47] is 
that there C — 0, and B = B(fi) is a scale dependent singlet function, while we define all entries 
of M as constants without specifying their relation to the scale dependent current correlators. 
In the isospin basis, (\(uu — dd) /y/2), \ {uu + dd) /y/2), \ ss)), M reads, 

(AO \ 

M — A + 2B y/2{B + C) , (76) 
V y/2{B + C) A + B + D J 

so that we can identify, A = M^ . The trace of M then yields the condition, 

W + D = Ml + Ml~ 2M p 2 , (77) 
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Since <C and T w <C M u , the complex phase in e can safely be neglected. 
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and e obtained in Appendix |A] gives the constraint, 



n B + C 

tan2e = \/8- -. (78 

D — _D 



The final relation, 

m - Mir 



3(B + C)-±(C + D) 



g(C + D)\ (79) 



shows that in the SU(3) limit, C = D = 0, singlet diagrams associated with B would split 
Mj from M 2 = Ml in much the same way as triangle anomaly diagrams would split M 2 , from 
M 2 = M 2 . Taking into account that e < 1, we can approximate, 

Ml - M\ _ M 2 + M% 

B * - 2 e , D » MJ - - 2 e , (80) 

which is correct up to 0(e 2 ). Thus, in the limit of ideal mixing, B drives the splitting of M^ 
from Mp instead. In any case, the singlet contribution associated to B is very small compared 
to A and D. More relevant for this work is the 5*^7(3) breaking singlet parameter, 

Mj - Ml Ml - M 2 n0 
* " tan2e p \ 81 



which reduces to C = —B in the ideal mixing case, e = 0. Numerically, we have, 

A = (769 MeV) 2 , B = (105 MeV) 2 , 
C = ( 74 MeV) 2 , D = (660 MeV) 2 . 



(82) 



There is a second solution in which \B\, \C\, and \D\ are all comparable and where B w — C 
to ensure e < 1. It also has D < 0, although the strange quark mass is expected to give 
the dominant (positive) contribution. Therefore we discard this solution. We can also roughly 
estimate the singlet component contained in D by relating it to M K *o, 

Aingict ~ D - 2{M 2 K , Q - M 2 ) = (123 MeV) 2 , (83) 

which is of similar size as B and C. Thus, singlet contributions to vector meson masses are 
generally suppressed by more than an order of magnitude beyond the QCD suppression factors 
discussed in Section [6] indicating that the smallness of the known coefficient, C2 C 1, may indeed 
be a generic feature that persists at higher orders. Qualitatively the same results are obtained 
if linear mass relations are used in place of masses squared, except that singlet contributions 
appear even more suppressed. 
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